Heavy-Quark Diffusion and Hadronization in Quark-Gluon Plasma 



Min He, 1 Raincr J. Fries, 1,2 and Ralf Rapp 1 

1 Cyclotron Institute and Department of Physics and Astronomy, 
Texas A&M University, College Station, TX 77843, USA 
2 RIKEN/BNL Research Center, Brookhaven National Laboratory, Upton, NY 11973, USA 

(Dated: June 30, 2011) 

We calculate diffusion and hadronization of heavy quarks in high-energy heavy-ion collisions 
implementing the notion of a strongly coupled quark-gluon plasma in both micro- and macroscopic 
components. The diffusion process is simulated using relativistic Fokker-Planck dynamics for elastic 
scattering in a hydrodynamic background. The heavy-quark transport coefficients in the medium 
are obtained from non-perturbative T-matrix interactions which build up resonant correlations close 
to the transition temperature. The latter also form the basis for hadronization of heavy quarks into 
heavy-flavor mesons via recombination with light quarks from the medium. The pertinent resonance 
recombination satisfies energy conservation and provides an equilibrium mapping between quark 
and meson distributions. The recombination probability is derived from the resonant heavy-quark 
scattering rate. Consequently, recombination dominates at low transverse momentum (jot) and 
yields to fragmentation at high pr- Our approach thus emphasizes the role of resonance correlations 
in the diffusion and hadronization processes. We calculate the nuclear modification factor and 
elliptic flow of D- and B-mesons for Au-Au collisions at the Relativistic Heavy Ion Collider, and 
compare their decay-electron spectra to available data. We also find that a realistic description of 
the medium flow is essential for a quantitative interpretation of the data. 
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I. INTRODUCTION 

Experiments at the Relativistic Heavy Ion Collider 
(RHIC) and the Large Hadron Collider (LHC) have been 
searching for the deconfincd phase of nuclear matter and 
have begun to probe its properties [HQ. There are strong 
indications that this new form of matter behaves like a 
nearly perfect fluid with high opacity and low viscos- 
ity, referred to as strongly coupled Quark-Gluon Plasma 

(sqgp) i,g. One of the major experimental findings 

is a large azimuthal anisotropy, V2 , in transverse momen- 
tum (pr) spectra of hadrons in non-central collisions [f|. 
To account for this observation, hydrodynamic simula- 
tions require an early initialization time,implying a rapid 
thermalization of the bulk medium However, the 

microscopic origin of the rapid thermalization remains a 
matter of debate. 

In contrast to light partons making up the bulk of 
the medium, heavy quarks (charm and bottom), pro- 
duced in primordial hard collisions and acting as im- 
purities in the QGP, are not expected to fully equili- 
brate with the surrounding medium. Due to their large 
masses (tuq) a memory of their interaction history may 
be preserved, thus providing a more direct probe of the 
medium properties than bulk observables [9l-tl2j. The 
thermal relaxation time of heavy quarks has been ar- 
gued to be larger than that of light quarks by a factor 
of m Q /T ps 5 - 20 [O, 03 (T: typical temperature of 
the QGP). As they diffuse through the medium, heavy 
quarks interact with th e lig ht partons and their spectrum 
becomes quenched [HI EJ- Moreover, as they couple to 
the collective flow of the medium in non-central heavy-ion 



collisions, heavy quarks may develop substantial momen- 
tum anisotropics. These two effects are translated into 
equivalent behavior of heavy-flavor (HF) meson (D and 
B) spectra and V2, and further into the spectrum and 
of their decay electrons. The latter have been measured 
in Au+Au collisions at RHIC fl3i - [l5j ). exhibiting appre- 
ciable modifications over their baseline spectra fromp+p 
and g?+Au collisions. 

Model calculations based on radiative energy loss in 
perturbative QCD (pQCD), which could account for the 
observed jet- quenching in the light sector [lj|, predicted 
a much smaller quenching for heavy quarks and asso- 
ciated single-electron spectra [13] ■ The large HQ mass 
suppresses small-angle gluon radiation ( "dead cone" ef- 
fect [lH) and reduces the gluon formation time [l9j . 
hence mitigating radiative energy loss significantly. How- 
ever, elastic collisions of heavy quarks with light par- 
tons [9i-[TTI. I20I l2lT | have been argued to dominate over 
radiative scattering at low momentum, resulting in no- 
table quenching of the HQ spectrum. 

However, jet quenching only captures part of the 
physics potential of the HQ probe. Its diffusion prop- 
erties, which reach all the way to zero momentum, in- 
clude energy-gain processes which are, e.g., instrumental 
for the coupling to the collective flow of the medium. 
Several studies of HQ diffusion have been conducted 
in recent years using Fokker-Planck [hJ HH, l2ll - |27j and 
Boltzmann transport [28l - [30T ] approaches, mostly imple- 
menting elastic collisions as the microscopic dynamics. 
They differ not only in their treatment of the background 
medium, but also in the evaluation of (a) the transport 
coefficients emerging from the interactions between the 
heavy quarks and the medium, and (b) hadronization of 
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heavy quarks into HF mesons. Concerning item (a), most 
studies employ variants of the pQCD interaction [3l| . 
while a novel approach with heavy-light resonant inter- 
actions was introduced in Refs. [Hi [22|. The latter was 
found to be a factor of 3-4 more efficient in HQ thermal- 
ization than pQCD, and was subsequently corroborated 
by microscopic T-matrix calculations u sing input poten- 
tials from lattice QCD (1QCD) [H [H Ht- Concerning 
HQ hadronization, several studies focused on indepen- 
dent fragmentation [25T4271 130| , which is not reliable in 
the low and intermediate-py regimes. Here, light par- 
tons surrounding the heavy quark have a high phase- 
space density which renders coalescence a more plausi- 
ble hadronization mechanism [35 38]. In Refs. (22T |23|. 
heavy-light quark recombination has been incorporated 
utilizing an instantaneous coalescence model [39( which 
could still be problematic at low pr due to lack of energy 
conservation. A reliable treatment of the 1ow-pt regime 
is important since the total number of heavy quarks is ex- 
pected to be conserved through the hadronization tran- 
sition. If the D- or B-meson spectra are distorted at 
low pt, the spectra at higher px are necessarily affected 
thus modifying the Raa (and V2) of D and B-mcsons 
and their decay electrons. 

The purpose of the present work is to establish a realis- 
tic and quantitative framework for HQ probes within (a) 
a strongly coupled QGP background medium (modeled 
by hydrodynamics), (b) a non-perturbative scenario of 
clastic diffusion in the QGP simulated by Fokkcr-Planck- 
Langcvin dynamics, and (c) a hadronization scheme at 
the phase transition based on the same interaction as in 
(b), combining recombination and fragmentation consis- 
tent with the limiting cases of kinetic equilibrium and 
vacuum hadronization. Unlike pre vious studies utilizing 
weak-coupling diffusion [TJ l24l. |26|. [27l |30| we try to im- 
plement the HQ probe consistently within a framework 
of strong coupling between heavy and light quarks, both 
in the QGP and during hadronization. Our comprehen- 
sive framework is hence conceptually compatible with the 
notion of a strongly interacting QGP. 

The strategy in this work is as follows. For the HQ 
transport coefficient we employ a non-perturbative T- 
matrix calculation of heavy-light quark interactions [HI, 
[33( | . This calculation supports Feshbach resonances in 
the QGP in the color-singlet and anti-triplet channels, 
surviving as rather broad states up to ~ 1.5 T c . They 
are responsible for the enhancement of the transport co- 
efficient compared to pQCD scattering. With these co- 
efficients we perform Langevin simulations of HQ diffu- 
sion through an expanding medium which is described by 
ideal 2+1-dimensional hydrodynamics (using the AZHY- 
DRO code [5(| at RHIC energies). At the phase tran- 
sition, heavy quarks are hadronized through coalescence 
with light quarks of the medium using the Resonance Re- 
combination Model (RRM) [4l| implemented on a hyper- 
surface given by the hydrodynamic simulation. The coa- 
lescence probability is evaluated using the resonant scat- 
tering rate of the heavy quark with light (anti) quarks, 



supplemented by independent fragmentation. The RRM 
formalism is consistent with the heavy-light Feshbach 
resonance formation found in the T-matrix used for the 
transport coefficient. This stipulates the role played by 
the resonance correlations in our work. With an arti- 
ficially large transport coefficient, we check the equilib- 
rium limit of the HQ distribution emanating from the 
combined hydro+Langevin simulation and the ensuing 
degree of equilibration of the HF mesons upon resonance 
recombination. The full space-momentum correlations 
generated by the hydro-Langevin simulation enter into 
resonance recombination. This enables a quantitative as- 
sessment of the radial medium flow on HF meson spectra 
at low pt as imprinted on the final Raa measurement. 

Our article is organized as follows. In Sec. |H] we in- 
troduce the ingredients for the hydro-Langevin simula- 
tion of HQ diffusion in the medium, i.e., the transport 
coefficient, the initial distribution in coordinate and mo- 
mentum space, and the background medium described 
by an ideal hydrodynamic model. Numerical results for 
the HQ Raa and V2 are discussed in the equilibrium limit 
as well as for realistic coefficients. Sec. Illll is devoted to 
HQ hadronization. We implement the RRM formalism 
on arbitrary hadronization hypersurfaces, elaborate the 
equilibrium mapping in resonance recombination, and de- 
termine the partition of coalescence and fragmentation. 
In Sec. IIVI we examine consequences of modifying the 
medium flow for the predicted HF meson spectra, trig- 
gered by indications that the partonic flow of the hydro- 
dynamic evolution is too soft. In Sec. [V] we make contact 
with current experiments in terms of the nuclear modi- 
fication factor and the elliptic flow of electrons from HF 
decays. In Sec. IVII we summarize and conclude. 

II. LANGEVIN SIMULATION OF HEAVY 
QUARK DIFFUSION 

A. Relativistic Langevin Kinetics 

The thermal momentum of a heavy quark at tem- 
peratures characteristic for heavy-ion collisions at RHIC 
amounts to ptn ~ y/rngT, which is parametrically larger 
than the typical momentum transfer, q ~ T, in a sin- 
gle elastic collision with a light parton from the bulk 
medium. Therefore many collisions are needed to change 
the HQ momentum considerably [13, El • This forms the 
basis for approximating the HQ motion in the QGP by a 
succession of uncorrelated momentum kicks and leads to 
a Fokkcr-Planck approach realized stochastically by the 
Langevin equations [jj [H, Hi| 

rfx = gdi, (1) 

dp = -T(p) pdt + y / 2D(p + dp) dt p , (2) 

where x and p are the position and momentum vector of 
the heavy quark, and E(p) = (trig + p 2 ) 1 ^ 2 is its energy. 
In the following we employ the post-point discretization 
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scheme in which the equilibrium condition (the relativis- 
tic fluctuation-dissipation theorem) takes the simple form 



D(p) = r(p) E{p) T 



(3) 



with T(p) being the drag coefficient and D(p) the (diag- 
onal) diffusion coefficient. The standard Gaussian noise 
variable, p, is distributed according to 



w(p) 



1 



(2tt) 3 / 2 



-p /2 



(4) 



Neither the original Fokker-Planck equation [12J, |42j nor 
the Langevin equation is Lorentz covariant. We choose 
the momentum and position updates for our HQ test par- 
ticles to be at equidistant time steps dr in the lab frame. 
For a flowing medium, as in our context, the momentum 
updates are rather to be done in the fluid rest frame. 
The updated 4-momentum is boosted back to the lab 
frame with the fluid four-velocity u^{x) = -f(v)(l,v(x)). 
The aforementioned equilibrium condition must be sat- 
isfied in order for the long-time limit of the test particle 
distribution to converge to the equilibrium (Boltzmann- 
Juttner) distribution as defined by the underlying back- 
ground medium. Further details of our algorithm will be 
detailed in a forthcoming article |44| . 



B. Thermal Relaxation Rate of Heavy Quarks 

The transport coefficient most commonly calculated 
from an underlying microscopic interaction of the heavy 
quark with the bulk medium is the thermal relaxation 
rate A(p; T). It is related to the drag coefficient, T(p; T), 
in the post-point Langevin scheme, Eq. ([2]), through 



T{p-T) = A{p-T) 



1 dD{p-T) 
E(p) 3E(p) 



(5) 



Utilizing the equilibrium condition ((3]) one can argue that 
T(p) = A(p) + 0{T/mq) and neglect terms to higher 
order in the inverse HQ mass (relative to the medium 
temperature). 



We employ HQ relaxation rates from Refs. |32j, where 
in-medium T-matrices have been calculated for both 
heavy-light and quarkonium channels. The input po- 
tentials were constructed using a field-theoretic ansatz 
for a confining and a color-Coulomb interaction with pa- 
rameters fitted to color-average free energies computed 
in finite-temperature lattice QCD (1QCD) This 
approach treats heavy quarkonia and heavy-light inter- 
actions in the QGP on an equal footing, and in both 
bound-state and scattering regimes. One thus obtains 
mutual constraints by analyzing, e.g., Euclidean corre- 
lation functions and HQ susceptibilities which turn out 
to agree fairly well with thermal 1QCD "data" [33j. For 
heavy-light quark scattering, the (non-perturbative) re- 
summation in the T-matrix generates resonances close to 
the 2-particle threshold (commonly referred to as "Fesh- 
bach resonances" ) in the attractive color-singlet (meson) 
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FIG. 1: (Color online) Imaginary part of the in-medium on- 
shell T-Matrix for charm-light quark scattering as a function 
of center-of-mass energy in the color-singlet (upper panel) 
and anti-triplet (lower panel) channels, taken from the lattice- 
QCD based potential approach of Ref. The vacuum T- 

matrices have been downscaled by a factor of 0.025. 



and color-anti-triplct (diquark) channels up to tempera- 
tures of about 1.5 T c , sec Fig.[T]for charm quarks (similar 
rcsults are obtained in the bottom sector). The increas- 
ing strength of the T-matrices in the color-singlet and 
anti-triplet channels when approaching T c from above 
is indicative for "pre-hadronic" correlations leading to 
hadronization. But even at high temperatures a sub- 
stantial enhancement of the T-matrix over elastic pQCD 
amplitudes persists, in particular close to threshold. 
The rather large resonance widths are mostly generated 
through the self-energies of the light- and heavy-quark 
propagators in the T-matrix (evaluated self-consistently 
in the HQ sector). 

The T-matrices have been used to calculate thermal re- 
laxation rates of heavy quarks [32l . HH . Resonant rescat- 
tcring accelerates kinetic equilibration by up to a fac- 
tor of ~3-5 relative to leading order (LO) pQCD calcu- 
lations @, cf. upper panel of Fig. [5] With increasing 
HQ 3-momcntum the thermal phase space of comoving 
partons (suitable for forming a Feshbach resonance) de- 
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FIG. 2: (Color online) (a) Charm-quark relaxation rate as a 
function of three-momentum using (i) heavy-light quark T- 
matrices (with 1QCD internal energy [45| as potential) plus 
pQCD gluon scattering with a a = 0.4 (upper 3 curves), 
and (ii) pQCD scattering off anti-/quarks and gluons with 
a s — 0.4 (lower 3 curves), (b) Temperature dependence of 
the charm/bottom quark thermal relaxation rate (at vanish- 
ing momentum) used in our simulations. The results are taken 
from Ref. [31. 



creases and the relaxation rate approaches the pQCD re- 
sults. For high energies and in Born approximation the 
T-matrix results recover the LO pQCD scattering ampli- 
tudes (32I ] . The temperature dependence of the charm 
and bottom relaxation rates (at vanishing 3-momcntum) 
used in our simulations is displayed in the lower panel of 
Fig. [2] They have been extrapolated linearly from the 
transition temperature in the 1QCD calculations of the 
free energies T c =196MeV, to T c =165MeV implicit 
in the equation of state as used in AZHYDRO. 



C. The hydrodynamic background QGP medium 

Hydrodynamic simulations are widely applied to model 
the bulk evolution of the matter created in heavy-ion col- 
lisions at RHIC providing a good description of 



hadron spectra and their elliptic flow. Here we use a 
hydrodynamic simulation of the fireball to provide the 
background medium for HQ diffusion. It supplies the 
information on the space-time evolution of energy and 
entropy density, as well as temperature and fluid veloc- 
ity which are needed to calculate the transport coeffi- 
cients in the Langevin dynamics and on the hadroniza- 
tion hypersurfacc. We have employed the publicly avail- 
able ideal 2+1-dimensional AZHYDRO code licjl in our 
study. It assumes longitudinal boost invariance |46| and 
has been tuned to fit to bulk observables at kinetic freeze- 
out at an energy density of ef = 0.075 GeV/fm 3 in 
v /^7 = 200 GeV Au+Au collisions at RHIC 0- The 
initialization of AZHYDRO is done at tq = 0.6 fm/c by 
specifying the entropy density distribution as 

1 3 

s(t , x, y; b) = K[-n B c(x,y, b) + -n WN (x, y; b)} , (6) 

where nee an d tiwn are the binary-collision and 
wounded-nucleon densities, respectively, calculated in the 
optical Glauber model [40| . and b is the impact param- 
eter. The coefficient n is fitted to the observed rapidity 
density of charged hadrons, dN^/dy, and translates into 
an initial entropy density of s(tq,0, 0;0) = 110/fm 3 at 
the center of the transverse plane for central Au+Au col- 
lisions at RHIC. 

In Fig. [3] we summarize the main features of AZHY- 
DRO relevant to our HQ diffusion calculations. The 
upper panel displays the time evolution of the energy- 
momentum anisotropy, e p = (T xx - T yy ) / (T xx + T yy ) 
for semi-central collisions (6=7 fm); it exhibits the de- 
velopment of the bulk anisotropy which leads to an el- 
liptic flow for final-state particles [iol] . One sees that 
e p tends to saturate at later times when the spatial 
anisotropy of the system has essentially vanished; the 
dip around r ~ 5 fm/c is due to the vanishing accel- 
eration in the mixed phase, which, in turn, is a result 
of the equation of state (EoS) with a Maxwell construc- 
tion between a non-interacting QGP with a bag constant, 
B = 0.3642 GcV/fm 3 at T > T c = 165 MeV, and a 
hadronic resonance gas at T < T c . Since in our Langevin 
simulations the HQ test particles freeze out at the end 
of the mixed phase (at edec = 0.445 GcV/fm 3 ), we show 
in the middle and lower panels of Fig. [3] the light-quark 
Pt-spectrum and V2 at this point, respectively. The light- 
quark mass is taken as m q = 350 MeV and we used 
the standard Cooper- Frye freeze-out procedure \4$, S3 ■ 
For comparison we also show the results of an empirical 
fireball parametrization of quark distributions extracted 
from multi-strange hadron spectra in Ref. [48| . The 
quark-p t spectra arc noticeably harder than in the hy- 
drodynamic evolution. Since multi-strange particles are 
believed to kinetically decouple close to T c , this suggests 
that the hydrodynamic evolution in the default AZHY- 
DRO docs not generate enough flow in the QGP. To in- 
vestigate the effect of a larger flow on HQ spectra we 
will also conduct Langevin simulations with a schematic 
fireball whose final-state flow is given by the empirically 
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FIG. 3: (Color online) (a) The time evolution of the asymme- 
try e p of the energy-momentum tensor in AZHYDRO for b = 7 
fm Au+Au collisions at t/snn — 200 GeV. (b) Light quark 
(m q = 350 MeV) pt-spectrum calculated with freeze-out at 
the end of the mixed phase with decoupling energy density 
e dcc = 0.445 GeV/fm 3 in AZHYDRO (red solid line). It is 
compared to the light quark spectrum at the end of the mixed 
phase of the parameterized elliptic fireball model discussed in 
Sec. IIVI (green dashed line), (c) The light quark elliptic flow 
V2 at the end of the mixed phase. Again, AZHYDRO and the 
parameterized fireball results are shown. 



extracted quark spectra. The pertinent elliptic flow ex- 
hibits slightly flatter p t dependence than in the hydro- 
dynamic simulation, cf. lower panel of Fig. G3 However, 
the integrated quark elliptic flow of (V2) = 4.99% is very 
close to the hydro result of (1)2) = 5.03%, representing 
the benchmarks from which the heavy quarks acquire V2 
through heavy-light parton interactions. However, an- 
other 20-30% is typically built up in the hadronic evolu- 
tion below T c (recall the upper panel of Fig. [3]) which is 
neglected in the present study. 



D. Initial Distributions of Heavy Quarks 

The number of heavy quarks produced in heavy-ion 
collisions is consistent with binary nucleon-nucleon col- 
lision scaling (49) . Thus their initial spatial distribu- 
tion is expected to follow the binary collision density, 
nBct^j U] b) which we adopt in our simulations within the 
transverse area where the energy density, e(ro,x, y), is 
larger than the decoupling value e^ cc — 0.445 GeV/fm 3 . 
For the initial HQ momentum distribution, we use the 
same spectrum as in Refs. [THHH, where PYTHIA re- 
sults for charm- and bottom-quark spectra, converted 
into D and B mesons via (5-function fragmentation, were 
tuned to semi-leptonic electron-decay spectra as mea- 
sured in p + p and (i+Au collisions at RHIC. This pro- 
cedure leads to a bottom-to-charm cross section ratio of 
a bll a cc = 4.9x 10~ 3 , and a crossing of the electron spec- 
tra from D- and i?-meson decays at pi ~ 5 GeV, see 
Fig. [U The b/c cross-section ratio is within the range of 
pQCD predictions [50l | and turns out to reproduce fairly 
well experimental data [HJ [HJ for the pt-dependence of 
the ratio of electrons from B-mesons to the sum from 
D+B, cf. lower panel of Fig. @] The i?-mcson contribu- 
tion becomes sizable for p\ > 3 GeV. 



E. Heavy-Quark Spectra and Elliptic Flow 

We now combine the ingredients as specified in the 
previous sections to perform the hydro+Langcvin simu- 
lation of HQ diffusion in the QGP using the test-particle 
method. A vector (xn,po) in transverse phase space, 
representing a heavy quark, is generated by Monte Carlo 
methods following the initial distributions discussed in 
Sec. Ill Dl Then we follow the trajectory of the heavy 
quark in phase space in equal time steps in the lab frame. 
At each time step, we read off the temperature, energy 
density and velocity of the fluid cell at the current HQ po- 
sition, (r, x, y, r] = 0). The drag coefficient is determined 
by the HQ momentum in the fluid rest frame and the tem- 
perature of the fluid cell. The momentum of the heavy 
quark is updated stochastically in the fluid rest frame ac- 
cording to the Langevin rule in Eq. ([2]) and boosted back 
to the lab frame using the fluid velocity. The HQ posi- 
tion is updated in the lab frame, which can be shown to 
be equivalent to an update in the fluid rest frame. Test 
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FIG. 4: (Color online) (a) Electron spectra from semileptonic 
decays of D- and _B-mesons (obtained from initial c- and b- 
quark spectra with 5-function fragmentation) in p+p collisions 
at RHIC energies, (b) Transverse-momentum dependence of 
the relative contribution of electrons from B-mesons to elec- 
trons from D+B decays. The solid curve results from the 
spectra in the upper panel which we adopt in our calculations; 
the data are from PHENIX [51| (filled squares) and from 
STAR [111 (filled circles) for p+p collisions at ^ = 200 GeV. 



FIG. 5: (Color online) (a) The charm quark pt-spectrum ob- 
tained from hydro+Langevin simulations with a large drag 
coefficient, F — 40.0/\/i?/fm (dots), compared to the equi- 
librated charm-quark spectrum calculated from the edcc = 
0.445 GeV/fm 3 freeze-out hypersurface in AZHYDRO (red 
solid line). The blue dashed line is the initial charm-quark 
spectrum with the same total yield, (b) The same compari- 
son as in (a) but for the elliptic flow. 



particles that have diffused away from r\ = to rapid- 
ity y and space-time rapidity r\ are redefined from the 
longitudinal phase-space coordinate (77; y) to (0; y — rf) to 
enforce boost-invariance. 

The heavy quark continues to diffuse in the QGP until 
the local energy density of the fluid drops below the de- 
coupling value, edcc = 0.445 GeV/fm 3 , corresponding to 
the end of the mixed phase of the cell. At that point, we 
assume the heavy quark to decouple from the fireball and 
mark it for hadronization. Wc do not take into account 
a possible local reheating if the expanding QGP phase 
"swallows" again an already hadronized heavy quark due 
to the increasing matter flow. Our criterion for the de- 
coupling of heavy quarks automatically yields their flux 
across the hadronization hypersurface as 

fQ(r,x,y;p)p IJ ,da fi (r,x,y)/E(p) (7) 



for any area element da^{T,x,y) on that surface, in ac- 
cordance with the Cooper-Frye formalism for the hydro- 
dynamic freeze-out. 

It is critical to verify that heavy quarks can reach 
local equilibrium as the stationary solution (53|. Wc 
have checked the equilibrium limit with an artificially 
increased drag coefficient, T = 40/ \JE [GcV]/fm, and 
a homogeneous initial spatial distribution for test parti- 
cles in the transverse plane. This specific choice for the 
energy dependence of T resembles the momentum de- 
pendence of the T- matrix based coefficients. 1 The size 
of the numerical coefficient (~40) in the large-T case is 



1 We have verified that the higher-order terms in Eq. (0, which 
are dropped in our Langevin simulations below, are negligible for 
the much smaller "realistic" coefficients. 
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limited by the requirement that the numerical time-step 
in the Langevin process be smaller than the inverse re- 
laxation rate. In the upper panel of Fig. [5] the Langevin 
charm-quark spectrum with large coefficients is compared 
to the distribution from Cooper-Frye freeze-out on the 
e doc = 0.445 GeV/fm 3 hypersurface in AZHYDRO, i.e., 
charm quarks in complete local thermal equilibrium. We 
have adopted a charm-quark mass of m c = 1.8 GeV, cor- 
responding to the in-medium mass at T c = 165 MeV 
in our simulations [HJ . The spectra agree well up to 
p t ~ 4.0 — 4.5 GeV. The deviation at higher p t is due 
to surface emission of charm quarks with large veloci- 
ties which escape the active (i.e. e(T,x,y) > edcc) part 
of the fireball at the earliest times; roughly 1% of heavy 
quarks at a given high p t do not suffer collisions, corre- 
sponding to the factor ~+00 suppression of spectra from 
the Langevin simulation relative to the initial distribu- 
tion at large pt- A matching picture is observed for the 
elliptic flow (lower panel in Fig. [5]): at low p t the v-i of 
the hydro+large-r-Langevin simulation follows the of 
equilibrated charm quarks, while it breaks away and os- 
cillates around zero for large pt (deviations set in slightly 
earlier than for the inclusive pt spectra, presumably since 
i)2 is a more differential and thus more "fragile" quan- 
tity). 

Next we turn to the results of our simulations un- 
der "realistic" conditions, using the transport coefficients 
and initial distributions outlined above, together with 
temperature-dependent in- medium HQ masses [32(. As 
usual, the modifications of the HQ spectra in the medium 
are quantified by the nuclear modification factor and el- 
liptic flow, 



RAA{Pt,y) 



V2(pt,y) 



dN AA 
dp t dy 



N co n 
I 



dN„ 



dptd(pdy COs(20) 



dN AA 



(8) 
(9) 



respectively, where iV co ii is the estimated number of bi- 
nary nucleon-nucleon collisions for the centrality bin un- 
der consideration. In Fig. [5] we display the charm- and 
bottom-quark Raa and vi at the end of the mixed 
phase as obtained from hydro+Langevin simulations 
in semi-central Au+Au collisions (6=7 fm). The ap- 
proach toward thcrmalization induces a depletion of 
heavy quarks at large pt (quenching) and an enhance- 
ment at low p t enforced by HQ number conservation. At 
Pt ~ 5 GeV, the charm-quark quenching reaches down 
to ~^0.4 while bottom quarks are much less affected, 
with i?AA(p*=5 GeV)~ 0.8. Note that at the same pt 
the Lorentz-7 of bottom quarks is significantly smaller 
than for charm. Radiative contributions to HQ transport 
are estimated to become competitive with elastic scat- 
tering once the non-perturbative effects are supp ressed, 
i.e., above p t — 4-5 GeV for charm quarks [12| (recall 
Fig. [5]). When using a drag coefficient from pQCD clastic 
scattering only (including both quarks and gluons with 
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FIG. 6: (Color online) Nuclear modification factor (upper 
panel) and elliptic flow (lower panel) of charm (red solid line) 
and bottom quarks (green dashed line) at hadronization ob- 
tained from hydro+Langevin simulations for 6=7 fm Au+Au 
collisions at RHIC energy, using transport coefficients from 
the heavy-light quark T-matrix plus a pQCD HQ-gluon con- 
tribution. For comparison charm-quark results are shown 
with coefficients using only LO pQCD scattering off gluons 
light quarks (blue dotted line). 



a s = 0.4), the quenching is weaker by about a factor of 
^3 [22[ • Similar features are found in the elliptic flow co- 
efficient, which for c-quarks first increases approximately 
linearly before leveling off at about 4.5%, characterizing 
a transition from a quasi-thermal to a kinetic regime. 

Previous calculations employing a thermal-fireball 
model for the medium evolution, using drag coefficients 
for non-perturbative elastic scattering [lol ] of compara- 
ble magnitude as in our calculation, have found signifi- 
cantly larger values for the maximal charm-quark V2 of 
around 7.5% Part of this difference originates from 
the larger "intrinsic" vi in the fireball medium which 
has been adjusted to the empirically observed hadron- 
V2 of 5.5-6%. Since the diffusion coefficient of charm 
(.D-mesons) in the hadronic phase is not negligible [H3 ] . 
the HF i>2 in the present study should be considered as 
a lower bound. Another source of uncertainty derives 
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from the freeze-out prescription and the associated real- 
ization of the HQ Langevin process (Cooper-Frye in the 
hydro evolution vs. Milckhin-likc in some fireball calcu- 
lations) [55| . 



III. HEAVY-QUARK HADRONIZATION 

The bulk matter in a hydrodynamic simulation can 
be evolved through a phase transition (here QGP to 
hadronic matter) solely by specifying the equation of 
state of the medium. However, the HQ spectra result- 
ing from the Langevin simulations through the QGP are, 
in general, not in full equilibrium with the bulk medium 
and thus require a microscopic hadronization mechanism 
to enable the calculation of HF observables. We will carry 
this out at the end of the mixed phase, represented by 
the hypersurface defined by the critical energy density of 
the hadronic phase in the hydrodynamic simulation. For 
simplicity we focus on the formation of D- and B-mesons 
neglecting HF baryons and hidden heavy flavor (both of 
which have been found to give small contributions to the 
total HF content of the hadronic phase [22| )■ Two mi- 
croscopic hadronization mechanisms have been consid- 
ered in heavy-ion physics to date: independent fragmen- 
tation of partons and coalescence of quarks. The former 
is appropriate for large-momentum partons emerging di- 
rectly from initial hard processes, with phenomenological 
fragmentation functions simulating vacuum gluon radia- 
tion and color neutralization. Coalescence, on the other 
hand, is believed to dominate in the low-momentum 
regime where partons are abundant in phase-space in 
heavy-ion I35 |-l37l and even in elementary hadronic re- 
actions j34 l56l|. 

Several previous studies of HQ diffusion in heavy- 
ion collisions have neglected coalescence processes |25l - 
l27l I30T ] , thus limiting the applicability of HF observables 
to high momenta. The formation time of heavy quarks 
is comparatively short, and thereafter their virtuality is 
small, governed by interactions with the medium with 
modest momentum transfers. Hence, fragmentation is 
not effective. In the Langevin simulations of Refs. [22I — 
|2~H , heavy- light quark recombination has been accounted 
for [3^] and found to be important for increasing both the 
elliptic flow and the nuclear modification factor of the 
resulting Z?-meson spectra. The coalescence formalism 
was based on the widely used instantaneous approxima- 
tion [35l [36j which, however, docs not conserve energy 
in the 2 — > 1 hadron formation process. A related prob- 
lem is the lack of a well-defined equilibrium limit for the 
hadron distributions. Together, both features imply ap- 
preciable uncertainties in calculating HF observables in 
the 1ow-pt region (albeit suppressed compared to light- 
quark coalescence by a mass ratio (m q /m c )). To improve 
the coalescence description and achieve consistency with 
local kinetic equilibrium, we here employ the resonance 
recombination model (RRM) implemented on the hydro- 
dynamic hadronization surface. 



A. Resonance Recombination at the Hadronization 
Hypersurface 

In the RRM the hadronization of constituent quarks 
is treated via resonance scattering within a Boltzmann 
transport equation [4l|. For scattering rates which 
are large compared to the inverse hadronization time, 
r rcs 1/Vhad, equilibrium quark distribution func- 
tions in a flowing medium are converted into equilib- 
rium meson spectra with the same flow properties, in- 
cluding elliptic anisotropics with space-momentum corre- 
lations characteristic for a hydrodynamically expanding 
source [HI, H3] • The RRM has been employed previously 
to investigate kinetic-energy and constituent-quark num- 
ber scaling [57} , and to extract empirical quark distribu- 
tion functions of the bulk medium at hadronization at 
RHIC (Hi. 

The RRM is consistent with the heavy-light Feshbach 
resonance formation found in the T-matrix calculation 
of the HQ thermal relaxation rate (see Section III Bp . It 
reiterates the important role played by resonance corre- 
lations in our work. As the temperature drops towards 
T c , the resonance correlations in the heavy-light quark 
T-matrix strengthen (recall Fig. [1]) and thus naturally 
merge into heavy-light quark recombination processes. 
When implementing the latter via a Breit-Wigner ansatz 
one obtains the HF meson distribution from the asymp- 
totic solution of the Boltzmann equation as [4l[ 



/M ymP (*,p) 



E M (p) 
m M T 



M 



d 3 p 1 d 3 p 2 , , . 
-^^/Q(x,Pi) 



X /g(x,p 2 ) O-(s) t>rel(Pl,P2) <5 (3) (P ~ Pi ~ P2) » (10) 

where /q, 9 ,m are equal-time phase-space distributions of 
heavy quarks, light quarks, mesons, respectively, 7J rc i is 
the relative velocity of the recombining heavy and light 
quarks, and 7Bm and Tm are the mass and width of 



the meson resonance [4lJ, |48|, |57|. In the calculations 
below we employ masses and widths compatible with 
the T-matrix calculation of Ref. (32j, extrapolated to 
T c = 165 MeV with m c = 1.8 GeV, m q = 0.35 GeV, 
m D = 2.25 GeV and T D = 0.1 GeV. 

Energy conservation and detailed balance in RRM en- 
sure an equilibrium mapping between the distributions 
of quarks and formed mesons [H, [5?} • We verify this 
in the present case of a non-trivial freeze-out hypersur- 
face given by AZHYDRO. We use local charm- and light- 
quark equilibrium phase-space distributions, f(p,x) = 
e - P -u(x)/T^ with fluid ve i 0C i t i es g i ven by AZHYDRO at 

the end of the mixed phase, then apply resonance re- 
combination, Eq. (jTUj) . locally (for each cell) to obtain 
the local meson phase-space distribution /m(t, x, y; p). 
Finally we calculate the current across the hypersurface 
and sum over all fluid cells on the edoc = 0.445 GeV/fm 3 
freeze-out hypersurface, 



dN 



PTdpxdtpdy 



P^d(J^(T,x,y) 
(2tt) 3 



/m(t,x,2/;p) . (11) 
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FIG. 7: (Color online) (a) D-meson pr-spectrum calculated 
with RRM on the AZHYDRO hadronization hypersurface 
(circles), compared to a direct calculation from AZHYDRO 
using the Cooper-Frye formula on the same hypersurface 
(solid line). The D-meson spectra at different rapidities 
(yD = 0.0 and yo = ±0.5) calculated from RRM agree with 
each other, (b) The same comparison for the elliptic flow of 
D mesons. 



In Fig. [7] we compare the resulting D-meson spectrum 
and V2 with a calculation directly from hydro using D- 
mcson Cooper-Fry freeze-out on the same hypersurface 
(edoc = 0.445 GeV/fm 3 ). The close agreement of the 
two calculations verifies the mapping between the equi- 
librium quark and meson distributions in RRM, includ- 
ing the full space-momentum correlations encoded in the 
AZHYDRO flow field. Longitudinal boost invariance of 
AZHYDRO is preserved by RRM as well, as observed 
from the independence of the D-mcson spectra on rapid- 
ity within our accuracy. 

The next step is to extend our approach to hadronize 
off-equilibrium quark distribution functions emerging 
from our HQ Langevin simulations. In order to couple 
the RRM to a HQ test particle freezing out from the 
hydro-Langevin simulation with momentum pdoc and co- 
ordinate Xdoc, we represent the corresponding local equal- 
time HQ phase-space distribution on the hadronization 



hypersurface, /q from Eq. 0, by a 5- function, 5 3 (x — 
Xdec)<5 3 (p- Pdoc) at the hadronization time T(x dc c, Udec)- 
As before, the light-quark phase space density f q at this 
point is taken to be the equilibrium distribution e~ p ' u / T 
at the local temperature and flow, and we can apply 
Eq. (flC)]) to obtain the phase space density /m of heavy 
mesons test particle by test particle. Finally, the spec- 
trum of heavy mesons follows from Eq. (| 1 1 [) as a sum 
over test particles. To check our procedure we first apply 
it to the Langevin output in the large-coefficient limit 
for charm quarks which - as discussed in Sec. Ill El - 
follows the equilibrium distribution up to pt — 4 GcV. 
Figure [8] shows the comparison between the D-meson 
spectrum and V2 calculated in this way and the D- 
mesons from a direct AZHYDRO calculation. Compared 
to Fig. [5j the agreement of the hydro+Langevin+RRM 
calculation with the direct hydro calculation extends to 
slightly larger px since the recombination essentially acts 
as an additional heavy-light interaction when forming D- 
mesons. This connection is particularly transparent if the 
interaction used in the diffusion process is the same as in 
the recombination process, i.e., a resonance interaction 
in our approach. We conclude that the equilibrium limit 
is well under control in this framework. 



B. Hadronization: Coalescence vs Fragmentation 

After establishing a coalescence formalism for 
hadronization of off-equilibrium HQ distributions on an 
arbitrary hypersurface it remains to couple this con- 
tribution with the standard fragmentation mechanism 
representing the largc-p t limit where the phase-space 
density of light partons vanishes (vacuum limit). In most 
previous works the coalescence probability has been 
evaluated in an instantaneous approximation which did 
not allow for a full control over its absolute magnitude. 
Here, instead, we use a dynamic criterion which directly 
follows from the RRM formalism; it is based on the 
HQ scattering rate which is derived from the same 
interactions as used in the diffusion calculations 2 . We 
thus make the following ansatz for the HQ coalescence 
probability in the fluid rest frame: 

P coal (p) = Ar rcs r r Q e » , (12) 

which is Lorentz invariant. In Eq. (|12p. the scattering 
rate, Tq s = n q (a q Q v Te \), refers to the resonant part of 
the Qq cross section, <r q Q (or T-matrix), and thus rep- 
resents the rate for hadron formation (n q : light-quark 
density, v rc i: relative velocity). The time interval Ar res 
characterizes the window in the dynamic medium evolu- 
tion during which resonance states exist; typically, this 



2 See Refs. \'2'i |54| | for a recent discussion of the relation between 
the two in the HF context 
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FIG. 8: (Color online) (a) D-meson pT-spectrum (stars) cal- 
culated from RRM on the hadronization hypersurface applied 
to charm-quark spectra from hydro+Langevin simulations in 
the large-drag coefficient limit (corresponding to Fig. [5} . It is 
compared to the D-meson spectrum directly calculated from 
AZHYDRO on the same hypersurface. (b) Same as in panel 
(a) but for D-meson elliptic flow. 



corresponds to the duration of the hadronization transi- 
tion, i.e., the "mixed phase", or even longer depending 
on whether "pre-resonance" states can be formed above 
T c . Of course, if the product Ar re s Fq s exceeds one, 
Pcoai should be put to one, corresponding to the equi- 
librium limit (more accurately, one could apply an ex- 
ponential relaxation, but in view of the practical uncer- 
tainties in the values for Ar ros and Tq s (j>) this is cur- 
rently not warranted). We emphasize that this proce- 
dure provides an absolute normalization of the coales- 
cence contribution, consistent with the (unique) equilib- 
rium limit. Since the resonance formation rate naturally 
diminishes with increasing HQ momentum (the phase- 
space density of quarks from the thermal bath to match 
the resonance mass decreases) , one obtains an increasing 
fraction, P[ rag (p) = 1 — Pcoai (p), of heavy quarks under- 
going independent fragmentation, recovering the vacuum 
limit. 

In practice, for our calculations reported below, we 



evaluate Eq. (|12|) as follows. For the HQ scattering 
rate we employ a Breit-Wigner cross section which is 
consistent with our heavy-light T-matrix (cf. Sec. Ill B[) 
and approximately reproduces the color-singlet contribu- 
tion to the HQ thermal relaxation rate at T c ; the perti- 
nent meson width of rM=0.4GeV is larger than the one 
used in the RRM expression, Eq. (|10j) . but the result- 
ing meson spectra and elliptic flow are rather insensitive 
to this quantity J4J [57} (we neglect resonant diquark 
contributions since we only consider color-singlet scat- 
tering relevant for D- and P-meson formation). With 
this cross section we calculate the HQ scattering rate in 
the fluid rest frame. We typically find T^ es « 0.1 GeV for 
charm quarks at vanishing momentum (similar for bot- 
tom quarks), consistent with Refs. [HI, [33| (about half 
of the total HQ width of ~0.2GeV calculated in these 
works is due to the color-singlet part). 

The HQ scattering rate is then boosted to the lab frame 
at the end of Langevin simulation (mixed phase) and ex- 
pressed as a function of HQ transverse momentum. For 
simplicity, we have chosen to apply P coa i not test-particle- 
by-test-particle but averaged over the spatial dependence 
in the fireball, i.e., as a function of pt only. We have 
checked that the explicit inclusion of space-momentum 
correlations leads to very similar results for heavy-meson 
spectra and 1)%. For the macroscopic time duration of 
resonance formation in the medium evolution, we adopt 
a conservative estimate of Ar rcs ~ 2fm/c (in the lab 
frame), amounting to P CO ai(pt — > 0) — >• 1; this is shorter 
than the duration of the mixed phase in the hydro evo- 
lution, Ar m j x ~ 3-4fm/c (roughly corresponding to the 
dip structure in the upper panel of Fig. [3]), and thus con- 
stitutes a lower limit of the coalescence contribution. 

In Fig. [9] we show the coalescence probability, P CO sd{pt), 
for charm and bottom quarks in the lab frame, aver- 
aged over all test particles with a given p t in semi-central 
(6=7 fm) Au+Au collisions at RHIC. At a given p t , bot- 
tom quarks have smaller velocities than charm quarks 
and thus more easily find a light-quark partner to recom- 
bine with. Consequently, the 6-quark coalescence proba- 
bility drops slower compared to c quarks. 

Let us summarize our procedure for HQ hadroniza- 
tion: for each HQ test particle we determine a heavy- 
meson spectrum for both recombination and fragmen- 
tation components, by applying either RRM or a frag- 
mentation function Dq(z) = S(z — 1). We then sum 
over test particles using respective weights P CO a\{pt) and 
1 — P C oai (pt) j yielding the full coalcscence-|- fragmentation 
meson spectrum. Thus, the total HF meson spectrum 
can then be written as 



dN tot & \ 



dN 



pxdpTdcjydy pxdpxdcfydy pxdpTd(f>dy 



(13) 



Note that we determine the absolute normalization of the 
spectrum by requiring the conservation of HF quantum 
numbers. Since our fragmentation function automati- 
cally conserves HF number we only need to impose the 
normalization on the recombination contribution to the 
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FIG. 9: (Color online) Charm- and bottom-quark coales- 
cence probability as a function of lab- frame pt for semi-central 
Au+Au collisions at RHIC. 



spectrum. We recall that the formation of heavy quarko- 
nia and heavy baryons has been neglected. 



C. Numerical results for D and B mesons 

Our final results for the D- and B-mcson pr-spectra 
and t)2 in semi-central Au+Au collisions at RHIC in the 
AZHYDRO+Langevin+RRM approach are displayed in 
Figs. ITD1 and [Til respectively. The absolute norm of the 
spectra is arbitrary (we have divided by the number of 
test particles, 5 x 10 7 , used in our simulations). We 
also show the individual recombination and fragmenta- 
tion components in the sum of Eq. (p~3|) . The coales- 
cence component dominates over fragmentation up to 
Pt — 3(7) GeV for D (B) mesons. Below these values, co- 
alescence with light quarks from the hydrodynamic heat 
bath increases the V2 by up to 30-40% compared to the 
HQ V2 (represented by the fragmentation component), 
whereas at high px the heavy-meson spectra and V2 ap- 
proach the fragmentation values. 

Rescattering effects due to HQ diffusion in the QGP 
are also clearly exhibited by the nuclear modification fac- 
tor which we display in Fig. [T2]for D- and .B-mesons (as 
well as for c and b quarks) in semi-central Au+Au colli- 
sions at RHIC. We note a dip of the D-meson Raa at low 
Pt, which is a consequence of the standard mass effect 
in the collective flow of the D-mesons. Remarkably, the 
dip is not present in the charm-quark Raa, underlining 
the effect of the larger D-meson mass and, more impor- 
tantly, the extra momentum added through coalescence 
with light quarks. In a sense, coalescence acts as an ad- 
ditional interaction driving the D-meson spectrum closer 
to equilibrium. The same effect is most likely responsible 
for the slight suppression of the D-mcson spectrum be- 
low the c-quark spectrum at high pt- Again, coalescence 
acts as an additional interaction towards equilibration, 
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FIG. 10: (Color online) (a) The coalescence, fragmentation, 
and total D-meson pr-spectrum for semi-central (6=7 fm) 
Au+Au collisions at y / %N = 200 GeV, normalized to one 
test particle, (b) The coalescence, fragmentation, and total 
elliptic flow of D-mesons from the same calculation. 

which in this case (high pr) leads to suppression relative 
to the c-quark spectrum. This is, in fact, contrary to 
the naive belief that coalescence should always add mo- 
mentum to the formed-hadron spectrum. However, the 
approach toward equilibrium is dictated by the underly- 
ing coalescence model (here RRM) possessing the correct 
equilibrium limit. 

The low-pT .D-meson spectrum is rather sensitive to 
the collective flow in the system. This is not the case for 
the D-meson spectrum as it does not come close enough 
to equilibrium, since the thermal 6-quark relaxation times 
are relatively large compared to the system's lifetime. In 
the following section we quantify how varying the collec- 
tive flow impacts charm (D-mesons) spectra. 

IV. MEDIUM FLOW EFFECT ON D-MESON 
SPECTRA 

In this section, we scrutinize the medium-flow ef- 
fect on D-meson spectra by comparing our results from 
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FIG. 11: (Color online) The same as Fig. llOl but for B-mesons. 
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FIG. 12: (Color online) D- and B-meson nuclear modification 
factors for semi-central (6=7 fm) Au+Au collisions at RHIC. 
For comparison, the charm and bottom quark Raa are also 
shown. 



the AZHYDRO background with those from a more 
schematic fireball model. We recall that AZHYDRO em- 
ploys an EoS with a mixed phase of appreciable time du- 



ration over which the sound velocity and thus the acceler- 
ation vanish. This is incompatible with state-of-the-art 
lattice QCD computations [581 ] . Consequently, AZHY- 
DRO presumably underestimates the flow at the end of 
the mixed phase and around chemical freeze-out. Rather, 
it has been tuned to fit bulk observables at kinetic freeze- 
out, Tf o =100McV; including multi-strange baryons such 
as the f2 _ • This is in conflict with the general belief 
(and empirical evidence) that multi-strange hadrons (f2, 
S, (j>) decouple close to T c (due to their small hadronic 
rescattering cross sections) but with rather large radial 
flow [H.[48ll59|. thus corroborating our assertion of insuf- 
ficient flow in AZHYDRO close to T c . 

Instead of retuning AZHYDRO with an improved EoS 
and/or initial flow [60(, we here adopt a modest atti- 
tude by investigating flow effects on D-meson spectra 
with a parameterized elliptic fireball model. We have 
modified the fireball model introduced in Ref. [22[ such 
that the light-quark p t -spectrum and its integrated el- 
liptic flow calculated at the end of the medium evolu- 
tion (QGP + mixed phase) agree with the empirical ex- 
traction of Ref. [48( where RRM was applied to exper- 
imental multi-strange hadron spectra and V2- For con- 
sistency with the hydro framework we adopt Cooper- 
Frye freeze-out rather than the Milckhin-like freeze-out in 
Ref. (cf. Ref. [Hj]). The retuned fireball results are in- 
cluded in Fig. [3] in direct comparison to those calculated 
in AZHYDRO. The integrated bulk V2 extracted from 
multi-strange particles (4.99%) is close to that calculated 
in AZHYDRO (5.03%) at T c , whereas the extracted light- 
quark spectrum is much harder than in AZHYDRO, due 
to the larger flow. 

We have utilized the fireball evolution to perform 
charm-quark Langevin simulations with our "realistic" 
coefficients and to compute £>-meson spectra from RRM, 
summarized in Fig. 1131 As expected, the D-meson px- 
spectra, shown in the upper panel for a hadronization 
with only recombination, are significantly harder com- 
pared to using the AZHYDRO background medium. The 
nuclear modification factor (middle panel of Fig. [T3")) ex- 
hibits a more pronounced flow effect at low momenta. 
The "flow-bump" for heavy particles is shifted to larger 
momenta in the fireball compared to AZHYDRO, and the 
depletion toward px = is larger. A significant part of 
this effect originates again from the coalescence process, 
as can be inferred from comparing the change from the c- 
quark to D-meson Raa in the "hard" fireball model (mid- 
dle panel of Fig|T3J), relative to the "softer" AZHYDRO 
calculation (Fig. [T2"|) , in the intermediate pr-region. This 
is, of course, due to the harder light-quark spectrum par- 
ticipating in heavy-light recombination. Also note again 
that the coalescence D-meson spectrum drops below the 
c-quark spectrum at high px, which once more illustrates 
the role of resonance recombination as an additional in- 
teraction driving the D-meson spectrum toward equilib- 
rium. The total D-meson V2, shown in the lower panel of 
FiglEl receives a modest increase for px > 3 GeV, due 
to a larger c-quark coalescence probability for stronger 
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FIG. 13: (Color online) Comparison of D-meson spectra and 
V2 from AZHYDRO and an elliptic fireball with larger flow in 
semi-central (6=7 fm) Au+Au collisions at ^/snn=200 GeV. 
(a) pT-spectra from resonance recombination only (normal- 
ized to the same total yield), (b) Nuclear modification fac- 
tor for /J-mesons (coalescence + fragmentation) and charm 
quarks obtained with the fireball model, (c) Elliptic flow 
from coalescence, fragmentation and their weighted sum for 
the fireball, and the total v 2 from AZHYDRO. 



FIG. 14: (Color online) (a) Electron Raa produced from semi- 
lcptonic D- and B-meson decays for central and semi-central 
Au+Au collisions using hydro and fireball background media, 
together with data from PHENIX for central and minimum 
bias collisions [3] an d from STAR for central collisions [liil ]. 
(b) The same for electron «2. 



V. HEAVY MESON SEMI-LEPTONIC DECAY 
AND OBSERVABLES 



Thus far at RHIC measurements of open HF in Au- 
Au collisions mostly pertain to their semi-leptonic single- 
electron decays fla - fla ]. 3 The latter have been shown 
to preserve the information on the nuclear modification 
factor and elliptic now of their parent hadrons rather 
well [H, HH . In the following we treat semi-leptonic de- 
cays of D- and S-mesons as free quark decays, c(b) — > 
s(c) + e + u e , albeit with effective quark masses equal 
to their mesonic bound states to correctly account for 
phase space: m;,=5.28GeV, m c =1.87GeV, to s =0.5 GeV, 
m e =0. 0005 GeV and to„=0 (light quarks are treated as 
spectators). We assume an average inclusive electronic 



flow (comoving light partons have a higher phase-space 
density at larger momentum when the flow is larger). 



3 Only very recently have direct D-meson data been reported. 



14 



branching ratio of 11.5% and 10.4% for c and 6, respec- 
tively. We have verified that hadronic form factors have 
little effect on the electron energy spectrum in the parent- 
particle rest frame, as already discussed in Ref. [62| . 
Specifically, the decays are performed via a Monte-Carlo 
simulation in the HQ rest frame, with the 3-body phase 
space weighted by the decay matrix element calculated in 
low-energy V- A theory [6l,[63j]: (\M\ 2 ) oc (p s -p v ){p c -p e ) 
for charm quarks, and (|A4| 2 ) oc (p c ■ p e )(Pb • Pv) for 
bottom quarks. Subsequently, the electron momentum 
is boosted to the lab frame using the calculated heavy- 
meson spectrum. 

Our numerical results for the single-electron Kkk(j>X) 
and V2(p1) in Au+Au at RHIC, based on the D- and 
B-spectra from our hydro+Langevin+RRM (as well as 
fireball+Langevin+RRM) calculations, are compared to 
data [HIHH m Fig-Ull For simplicity, we adopted an im- 
pact parameter of 6=7(0) fm to mimic the experimental 
minimum-bias (central) results. On the one hand, the 
minimum bias sample for bulk observables is closer to 
6=8-8.5 fm, but, on the other hand, HQ production scales 
more strongly with centrality, ^A 4 / 3 , implying a smaller 
impact parameter for the HF minimum bias sample. Our 
"conservative" choice of somewhat smaller 6 could thus 
slightly overestimate the suppression and underestimate 
the V2- 

In the upper panel of Fig. [TJ] one sees that the hydro- 
based calculations overestimate the suppression found in 
the experimental data for R e AA in the regime where the 
uncertainty of the latter is relatively small. The fire- 
ball calculations with larger (and probably more real- 
istic) flow improve on this aspect, which reiterates the 
importance of medium collectivity in the recombination 
process. For the elliptic flow (lower panel of Fig. IT4"|) wc 
find good agreement of both calculations with PHENIX 
run- 7 data up to pi ~ 1.5 GeV, whereas the run-4 data 
are underpredicted. At larger p\ our calculations tend to 
underestimate the run-7 «2 data. This indicates that ad- 
ditional contributions to charm and bottom interactions 
are required, e.g., non-perturbative HQ-gluon scattering 
(presently we treat this part perturbatively), radiative 
scatte ring and diffusion of HF hadrons in the hadronic 
phase [l2|. We have recently estimated the charm diffu- 
sion coefficient in the dense hadronic phase to be compa- 
rable to the values in QGP that we use here [54 1. 



VI. SUMMARY AND CONCLUSION 

In this work, we have developed a framework in which 
heavy-quark diffusion and hadronization in a quark gluon 
plasma are evaluated consistently in a strongly coupled 
(non-perturbative) scenario. The strong coupling is re- 
alized by resonance correlations which build up in the 
hadronic channels of anti- / quark correlation functions as 
the system cools down toward T c . On the one hand, 
HQ transport has been based on heavy-light quark T- 
matriccs (consistent with vacuum spectroscopy and the 



pQCD limit of high-energy scattering), in which meson 
and diquark resonances enhance the HQ relaxation rate 
over perturbative calculations [H, [H| ; these coefficients 
have been utilized in relativistic Langevin simulations 
of HQ diffusion with a medium evolution described by 
an ideal hydrodynamic model (which itself is based on 
the strong-coupling limit). On the other hand, as the 
medium temperature drops towards the critical value, 
the resonance correlations strengthen and naturally trig- 
ger heavy-light quark recombination, which is carried out 
in the RRM formalism. We have verified the equilib- 
rium mapping between quark- and meson-distributions 
on non-trivial hadronization hypersurfaces in the hydro- 
dynamic medium, which, in particular, allowed us to 
identify and quantify the important effect of the medium 
flow on the D-meson spectrum. We have also given a 
more rigorous definition of the coalescence probability 
in terms of the underlying formation rate, and in this 
way determined the partition between recombination and 
fragmentation contributions to the heavy-light spectra in 
absolute terms. We recover the mandatory limits of equi- 
librium and independent fragmentation for low and large 
transverse momenta, respectively. 

Despite a few missing components in our HQ-transport 
and hadronization scheme (as discussed below), we have 
carried our calculations to the level of HF observabels in 
heavy-ion collisions. In particular, wc predict that the 
degree of charm-quark thcrmalization is large enough to 
develop a characteristic "flow-bump" in the D-meson nu- 
clear modification factor. The proper equilibrium limit 
of the underlying coalescence mechanism is essential in 
developing this feature, while its location in px is sensi- 
tive to the strength of the medium flow. For 6-quarks 
the coupling to the medium is not strong enough for 
this feature to emerge. We have further decayed our D- 
and i?-meson spectra semi-leptonically; the correspond- 
ing single-electron spectra show an encouraging agree- 
ment with current RHIC data for p\ <1 GeV; the nuclear 
modification factor reiterates the importance of a realis- 
tic flow strength, while the elliptic flow is a more sensitive 
gauge for the magnitude of the HQ transport coefficient. 
Our current estimate confirms earlier results for the spa- 
tial diffusion constant in the vicinity of 5/ (27rT) . 

Several uncertainties and areas of improvement re- 
main, e.g., (i) a more complete evaluation of HQ re- 
laxation rates by including non-perturbative effects in 
elastic scattering off gluons, as well as adding radiative 
processes (expected to become relevant at high p t ); (ii) 
hadronic diffusion, especially in light of recent estimates 
indicating comparable strength of HF transport in the 
hadronic and QGP phases near T c (&4| ; (iii) a more com- 
plete hadro-chemistry at T c including strange _D-mesons 
and charmed baryons; (iv) a possible Cronin effect in the 
initial HQ spectra, which mostly affects Raa', (v) further 
improvements in the medium evolution, i.e., a hydrody- 
namic background with harder spectra (e.g., due to vis- 
cosity, initial flow or a harder EoS). The framework put 
forth in this paper allows us to systematically address 
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all of these items. This will not only be beneficial for 
the open heavy-flavor sector (and an ultimate determi- 
nation of the QGP's viscosity) but also for its impact on 
closely related observables such as heavy quarkonia and 
intermediate-mass dileptons. 
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